Particle addition and subtraction as a test of bosonic quality 
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We propose a test to measure the bosonic quality of particles with respect to physical operations 
of single-particle addition and subtraction. We apply our test to investigate bosonic properties of 
composite particles made of an even number of fermions and suggest its experimental implementa- 
tion. Furthermore, we discuss the features of the processes of particle addition and subtraction in 
terms of optimal quantum operations. 
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Introduction. Non-commutativity of photonic addition 
and subtraction was recently confirmed in the laboratory 
[HQ- An initial thermal state of light was subjected to 
both single photon addition and subtraction, mathemat- 
ically denoted by bosonic creation and annihilation op- 
erators, a) and a, respectively. Two scenarios were con- 
sidered, photonic addition followed by subtraction (AS) 
and vice versa (SA). According to quantum theory, the 
final states in the two scenarios should differ due to the 
bosonic commutation rule [a, a'] = 1. In general 

Psa = {a)a)p(aa)) = V p n , m nm\n){m\ ^ 



(1) 



Pas = (aa^)p(a r a) 



tn ^ = Pn,m(n + l)(m + l)\n)(m\ 



where \n) denotes the state of n indistinguishable parti- 
cles in the same mode; for simplicity the states are not 
normalized. 

Here, we show that bosonic addition and subtraction 
can be used not only to verify the commutation rela- 
tions of ideal bosons, but also to detect and infer some 
structural properties of composite particles (consisting of 
an even number of fermions) as well as to test to what 
extent they behave like ideal bosons. The paper is or- 
ganized as follows. First, we introduce a simple test to 
measure bosonic quality based on single particle addi- 
tion followed by subtraction. Since we are not directly 
interested in testing commutation relations we simplify 
the previous scenarios and do not consider the process 
in which subtraction is followed by addition. This also 
allows us to bypass one important problem, namely that 
the operations of addition and subtraction do not com- 
mute even for classical distinguishable particles. We then 
discuss general quantum operations (channels) that ef- 
ficiently add and subtract a particle. We identify im- 
portant criteria that have to be satisfied by bosonic ad- 
dition/subtraction channels and apply these channels to 
composite particles. We find that our measure of bosonic 
quality depends on the entanglement between fermions, 
therefore it provides us with information about the inter- 
nal structure of the composite particle and can be used as 



an entanglement witness. We conclude with a discussion 
of how to implement bosonic addition and subtraction of 
composite particles in the laboratory. 

Measure of bosonic quality. In this work, we pro- 
pose to compare the resulting state after particle ad- 
dition followed by subtraction with the initial state of 
the particles. Note that the operation AS for distin- 
guishable particles, described by the creation operator 
a} d = X)^Lo l n + l)( n l an d the corresponding annihilation 
operator, leaves the initial state unchanged unlike in the 
case of ideal bosons. However, the action of SA changes 
even the state of distinguishable particles. A successful 
subtraction indicates the absence of vacuum in the initial 
state and one has to update and re-normalize the state of 
the system. This lower bound on the number of particles 
is the heart of the classical lack of commutation of par- 
ticle addition and subtraction [a^a^] = |0)(0|. In order 
to get rid of this classical component in a test of bosonic 
commutation relations, one has to consider states that do 
not contain any vacuum. Since both states in [l], Q were 
thermal states having a substantial occupancy of vacuum 
(exceeding 50%), it would be interesting to repeat these 
experiments on vacuum-less states. 

Apart from the consideration of just the AS operation 
and the initial state, we make one further simplification 
for the sake of experimental convenience. Note that in 
general, to compare two states (for example using the fi- 
delity measure as done in [H0]), it is necessary to perform 
state tomography which is experimentally an extremely 
demanding process. The calculation of fidelity is how- 
ever not necessary to detect the change caused by AS, 
it is sufficient to measure the probability distribution of 
the number of particles {pa,pi, ■ ■ .}, where p n denotes 
the probability of detecting n particles. This is still com- 
plicated due to the fact that the current particle detec- 
tors do not efficiently measure the number of particles. 
Instead, their action can be described by measurement 
operators {|0)(0|,1 — |0)(0|} (see additional material in 
PJ). These detectors allow the measurement of po, which 
is sufficient to detect the change due to AS, and conse- 
quently, to provide a measure of bosonic quality. 
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The action of bosonic AS affects the probability dis- 
tribution in the following manner: p n — > (n + l) 2 p n , 
which together with normalization implies a decrease in 
Po- In contrast, the AS operation for distinguishable par- 
ticles leaves po unchanged. Due to the normalization the 
change in po depends on the total probability distribution 
{p n }- Let us now examine the difference in po, 



A <? 

PO ~ Po = PO 



PO 



= P0 



Po 



((N + l) 2 ) ' 



where p^ s denotes the vacuum occupancy after AS, N 
denotes the particle number operator and n max denotes 
the maximum number of particles in the system. Simple 
optimization yields that the maximal change occurs for 
Po — ^ max ^2 and Pn max = 1 — Po- The greater the max- 
imal number of particles n max , the greater the change 
in po after AS. At this point, it is clear that the above 
difference can be taken as a measure of bosonic quality 
M = po — Po S ■ F° r distinguishable particles this is al- 
ways zero, whereas for true bosons this is = n nma + 2 » 
approaching one in the limit of large n max . We propose 
to consider the simplest optimal probability distribution 
taking into account only the vacuum and a single particle 
component {po = §>f>i = ■§}, for which Mb = |. Note 
that the mixture po|0}(0| + (1 — po)|l)(l| can be prepared 
for any type of particle (even for fermions). The values 
M = and M = | correspond to distinguishable par- 
ticles and bosons, respectively. For simplicity, wc define 
our measure of bosonic quality as 



M 



3p£ 



AS 



(2) 



For ideal bosons, Mb = 1 and for distinguishable parti- 
cles Aid — 0. Any value between and 1 corresponds to 
imperfect bosons. 

Addition and subtraction channels. We now move 
to some general considerations of the particle addition 
and subtraction operations. In general, the processes 
of particle addition and subtraction are not determin- 
istic. Moreover, they cannot be formulated simply as 
Kraus operators K j which describe nondeterministic evo- 
lutions in terms of completely positive quantum channels 
p' = J^j KjPK*j 3]. The reason is that all the Kraus op- 
erators {-Koj K\, . . . } which describe a quantum channel 
cannot increase the norm of the state, i.e., J^j KjKj < 1. 



N + l. Theeigen- 



Setting K = a 1 yields K^Ko = at 
values of the bosonic particle number operator N lie in 
the set of all nonnegative integers, which together with 
the requirement that the norm cannot increase immedi- 
ately implies the negativity of the remaining operators 
K]Kj for j 0. The case of the annihilation process is 
analogous. It is interesting to note that the probabilis- 
tic nature of the operators a' and a can also be deduced 
from the fact that deterministic addition and subtraction 
could lead to an increase of entanglement via local oper- 
ations. This can be seen for example by considering the 



state of a single particle in two modes A and B, = 
a|CUl,B) +/3\Ia0b) with real parameters a 2 + f3 2 = 1 
and a 2 > 2/3 2 , that has entanglement measured by con- 
currence given as 2a(3. It is clear that a local operation of 
addition followed by subtraction at mode A leads to the 
state \ipAs) = a ilipi (c*|0a1b) + 2/3|1aO b )) with entan- 

which 
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glement measured by concurrence given as — 
is larger than the initial entanglement. The probabilistic 
nature of the operators a> and a is needed to ensure that 
entanglement does not increase via local operations. 

Note that any operator which effectively adds one par- 
ticle to the system is of the form 



"eff 



n=0 



fn\n+l)(n\ 



(3) 



An effective annihilation operator is the hermitian conju- 
gate of the above. This operator is a valid Kraus operator 
if \fn\ 2 < 1 for all n. It is convenient to rewrite it in the 
following form 



g{N)a) = ^2g(n + l)Vn + l\n + l){n\, 

n=0 



where g(N) is a function of the particle number opera- 
tor. The extreme case in which all multiplicative factors 
are equal to one corresponds to the creation operator 
of distinguishable particles at, where g(N) = 1/yfN. 
However, the test requires an operation that changes 
the ratios between probabilities p n , which can only be 
done if the multiplicative factors are different for each 
n. Optimally, we would like to have the ratios to be 
= ^^-t 1 . In this case g(N) is a constant, how- 
ever the normalization constraint and the fact that the 
sum over n goes to infinity imply that the only pos- 
sible solution is trivial g(N) — 0. This problem can 
be circumvented if the maximal number of particles is 
bounded. In this case, the optimal operator a e ff is state 
dependent, i.e., for the state supported on the subspace 
spanned by {|0), |1), . . . , \ n max }} the corresponding func- 
tion is g(N) = ^ n = +1 , a constant for n < n max , and 

g(N) — -j= for n > n max . The effective operator 



for n > n„ u 

with this function is then the optimal operator to imple- 
ment particle addition. 

Composite particles of two distinguishable fermions. 
We now proceed to situations for which the measure M. 
lies between and 1. Apart from any experimental im- 
perfections, this may occur if one deals with composite 
particles, i.e., systems which are composed of an even 
number of elementary fermions. Here, we discuss the 
case of a composite particle made of two distinguishable 
elementary fermions (see Refs. [iMZj]), the simplest de- 
parture from ideal bosonic behavior. 

The general pure state of two distinguishable fermions 
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can be written as 



AB=J2^a{b{\0), 

k 



(4) 



where at (frj.) creates a fermion A (B) in mode k and 
Afe are probabilities. The above state is written in the 
Schmidt form, i.e., as a sum of tensor products of dis- 
tinct orthogonal states. The modes k can refer for in- 
stance to energy levels of a confining potential, or to the 
position of the center of mass of A and B. The state 
can be considered as a single boson state if the operator 



c* = \/Aka k b k behaves as a proper bosonic creation 
operator. The corresponding commutation relation reads 
[c, ct] = 1 — A, where A = J^k ^fc( a ! a fc The state 



of n composite bosons is 



\n) = X, 



-1/2 _£ 



t" 



|0>, 



on which the action of the annihilation operator gives 
c\n) = a n \/n\n — 1) + \£n)- The parameters \n and 
OL n = \/xn/Xn-i are normalization constants and |e„) is 
a vector orthogonal to |n — 1). Moreover, 

(e n \e n ) = 1 -„-*-_ + ( n -l)-2±i. ( 5 ) 

Xra-1 Xn 

In the ideal scenario a n — > 1 and (e n \e n ) — > 0, which 
happens only if the ratio Xn±i/Xn ~> 1- This ratio was 
related to the parameters A& and it was shown that it 
approaches one in the limit of infinite entanglement be- 
tween A and B 

Now, let us examine the action of particle addition 
for composite particles. One can always write the corre- 
sponding Kraus channel as 



E 

71 = 



g(n + l)a n +iVn+ l|n + l)(n|, 



where the terms related to |e„) are incorporated in other 
Kraus operators. The particle subtraction is given by 
taking hermitian conjugate of the above. The optimal 
function is again a constant, g(n + l) = — 



for which the ratio 



fn 



.„ a , + l\Aw+l' 

The action of 



addition and subtraction on the optimal initial state 



f|0)(0| + i|l)(l| gives 



deduce the ratio x ^ +1 from A4. Repeating this for all n 
one can obtain full information about the internal state 
of the composite particle. We note that a high value 
of measure A4 for a composite particle implies that the 
deviation from the bosonic commutation is minimized 
for states \n) i.e., (n\A\n) scales inversely as M.. The 
measure M is thus an indicator of the bosonic quality 
of composite particles not only with respect to addition 
and subtraction but also with respect to other bosonic 
properties such as condensation, bunching etc. 

It is interesting to test our measure on the example of 
maximally entangled states with Afe = 4, where d is the 
total number of modes in the system. The corresponding 
creation operator is 



k=l 



Note that the consideration of a finite dimension d is 
natural when one deals with excitations confined in a 
finite condensed matter region. Moreover, if one deals 
with particles in free space, due to the finiteness of phys- 
ical resources such as energy, it is natural to assume fi- 
nite sums over coefficients A& . In the set of all operators 
spanned over d fermionic modes, operators (y max are the 
most bosonic since they correspond to maximally entan- 
gled states. Therefore, the bosonic quality of all states in 
such a set is bounded from above by the bosonic quality of 
the state defined by c\ nax . For these truncated operators, 
one is always restricted to states \n) for which n < d. It 



|0) 



\(d-n)\ 



hence 



is easy to verify that (0|c™ ax c% ax 
Xn = d\/ (d n (d — n)\) and a n = (d — n+ l)/d, which im- 
plies that the norm in ([5]) is zero. The highest quality of 
bosons obtainable with composite particles is 



M = 1- 



2d- 1 



3d 2 -2d+l 



(6) 



In the limit of infinite d (infinite entanglement) the above 
measure goes to one, which agrees with the results pre- 
sented in the Refs. [B-Q]. 

Composite particles of an even number of fermions. 
A general composite particle made of an even number 
(2n) of fermions, such as an atom or an ion, has creation 
operator given by 



multi 



'Pjk e jk+i 



■ e ]i n ]i+i 



M 



1 



2x1 ■ 



Therefore, one can obtain \i from this test. The purity 
P = Tr (Tr b(pab)) 2 , which is a measure of entanglement 
for pure bipartite states pab, is given by P = 1 — X2, 
and therefore the above test can detect entanglement be- 
tween the two constituent fermions. For initial states of 
the form po|0)(0| + (1 — po)\n)(n\, one can immediately 



where operators p\ , and n k may be taken to correspond 
to protons, electrons and neutrons in modes i, j and k 
respectively and rjj is a complex coefficient for a partic- 
ular configuration j — {ji, . . . , j2n}- It should be noted 
that the Kraus channel for addition takes an analogous 
form to Eq. ((3]) and the value of the measure M. is still 
a good indicator of the bosonic quality of the composite 
particle. 
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As regards the entanglement structure of these com- 
posite particles, note that this entanglement is tripartite 
with respect to the groups of particles of types p, e and 
n. Nevertheless, the purity of the reduced density matrix 
of any single subsystem is a measure of the entanglement 
between this subsystem and the rest of the particles. For 
genuine tripartite entanglement of the GHZ type [1] , the 
measure M can still be used to calculate the entangle- 
ment in the state. It is an interesting question how much 
information about other kinds of tripartite entanglement 
is contained in M. 

Implementation. So far, we have considered addition 
and subtraction operations implemented by the optimal 
Kraus channels. We now discuss possible experimen- 
tal implementations of these channels. In experiments 
P, [H, photonic addition was implemented by the pro- 
cess of parametric down-conversion (PDC), while pho- 
tonic subtraction was realized using a beam splitter (BS) 
with very low reflectivity. In both instances, successful 
implementation was heralded by a detection event on an 
auxiliary detector — for PDC a trigger photon confirmed 
the addition of the idler photon to the state, whereas for 
BS the heralding was done by the reflected photon sub- 
tracted from the state. The effective addition operator 
in the process of PDC is 

oo 

a PDC — 51 s'm^Wn + l)\n + l)(n|, 

where 7 is the PDC efficiency and t is the time of inter- 
action between the PDC crystal and the pumping field 
photon. For small 7 and t, (t is of the order l/c, where I 
is the size of the PDC crystal and c is the speed of light 
in the crystal) the above process well approximates the 
optimal bosonic addition channel. Moreover, the effec- 
tive subtraction operator resulting from a beam splitter 
of low reflectivity (r <C 1) is of the form 

00 

a BS = ^/{l-r)r^/n+l\n) (n + l\, 

which approximates bosonic subtraction very well. 

For massive particles, a possible candidate to realize 
the optimal addition channel is the process of Feshbach 
resonance, whose effective Hamiltonian is of the form 
Hf = r y(b ji aia2+ba\al), where a \ creates an atom of type 
1, a\ an atom of type 2 and a two-atom molecule. The 
optimal subtraction channel for massive particles can be 
implemented with two potential wells for which the tun- 
neling probability from one to the other is small, effec- 
tively mirroring the beam splitter operation. Note that 
these operations are in principle implcmcntable with cur- 
rent technology with Cs and Rb atoms |9[ . 

Discussion. The measure M. depends on two crucial 
features of the addition and subtraction process, namely, 
the quality of the particle that is to be determined by 



the test and the quality of the channels used to imple- 
ment these operations (i.e., how well they approximate 
the optimal addition/subtraction channel). For the opti- 
mal channels, as discussed previously, a good initial state 
is po = I and p\ = | and the measure takes the form © . 
In general, the value of the measure changes with the f n 
in the channels that implement the AS operation. Since 
the test is aimed at determining the quality of particles 
irrespective of the channel used, we now discuss a way to 
remove the effects due to imperfect channels. 

For general channels with the fraction -J 2 ^- ^ ^ n + 1 it 
can be easily verified that the measure reads 




The optimal initial state to maximize the value of this 

_ e , ll e„_«h e f„™ w .l/(l + ^) 2 ) 

and pi = 1 — pq. Therefore, while the optimal initial state 
differs for different channels the measure M. is always 
a valid indicator of bosonic quality irrespective of the 
addition/subtraction channel used. 

Conclusions. In this paper we proposed a test to deter- 
mine the bosonic quality of particles by elementary op- 
erations of single particle addition and subtraction. We 
discussed the probabilistic nature of these operations by 
formulating them within the Kraus operator formalism. 
By comparing an initial state with one after the process 
of addition followed by subtraction we formulated a mea- 
sure that determines the bosonic quality of the particle 
independently of the quality of channels performing these 
operations. We also suggested an experimental realiza- 
tion of these channels for composite particles. 
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